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^ ■ ABSTRACT 

^ . Local scaling properties of the co-added foreground-cleaned three-year Wilkinson 

^T) ' Microwave Anisotropy Probe (WMAP) data are estimated using weighted scaling in- 

\^ [ dices a. The scaling index method (SIM) is - for the first time - adapted and applied 

to the case of spherical symmetric spatial data. The results are compared with 1000 
Cn| ' Monte Carlo simulations based on Gaussian fluctuations with a best fit ACDM power 

, spectrum and WMAP-like beam and noise properties. Statistical quantities based on 

the scaling indices, namely the moments of the distribution and probability-based 
measures are determined. We find for most of the test statistics significant deviations 
from the Gaussian hypothesis. Using a very conservative statistics, which aver- 
pi i' ages over all scales, we detect non-Gaussianity with a probability of 97.3 % regarding 

the KpO-masked full sky, 98.9 % for the KpO-masked northern and 91.6 % for the 
KpO-masked southern hemisphere. When analysing different length scales separately, 
the detection rates increase to 99.7 % (> 99.9 % north, 97.4 % south) for the mean, 
^ '. 98.4 % (99.9 % north, 71.6 % south) for the standard deviation and 97.9 % (99.4 

% north, 80.3 % south) for a x^-combination of mean and standard deviation. We 
find pronounced asymmetries, which can be interpreted as a global lack of structure 
in the northern hemisphere, which is consistent with previous findings. Furthermore, 
''^Jj ' we detect a localized anomaly in the southern hemisphere, which gives rise to highly 

I significant signature for non-Gaussianity in the spectrum of scaling indices P{a). We 

identify this signature as the cold spot, which was also already detected in the first 
year WMAP data. Our results provide further evidence for both the presence of non- 
Gaussianities and asymmetries in the WMAP three-year data. More detailed band- 
and year- wise analyses are needed to elucidate the origin of the detected anomalies. 
In either case the scaling indices provide powerful nonlinear statistics to analyse CMB 
maps. 

Key words: cosmic microwave background - cosmology: observations - methods: 
data analysis 
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1 INTRODUCTION 

The study of the Gaussianity of the Cosmic Microwave Back- 
ground (CMB) is regarded to be the best way for under- 
standing the true nature of the primordial density fluctu- 
ations; by measuring suitable statistics of the temperature 
fluctuations of the CMB and comparing the results with the- 
oretical predictions, cosmological models for the primordial 
density fluctuations can be rejected or accepted at a cer- 
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tain confidence level. Standard inflationary m odels (jCuthl 
Il98ll : lLindelll982l : lAlbrecht fc Steinhard|[l983 ) predict that 
the temperature fluctuations of the CMB correspond to a 
(nearly) Gaussian, homogeneous and isotropic random field. 
In fact, Gaussianity and statistical isotropy are among the 
fundamental pillars of the ACDM concordance cosmologi- 
cal model. On the other hand, many alternative scenarios 
have been studied, whi ch give rise to non-Gaussianity, e.g. 
non-standard inflation (iLinde fc Mukhanov I [l997l: iPeebled 



1 1997': 'Bcr nardeau fc Uzanll200j: lAcauaviva et al.H2003l). or 
topological defects models ( Bouchet. Bennett, fc Stebbind 
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1 19881 : iTurok. fc Spergel I Il990l : iTurokl Il996l ) to mention a 
few. Probing non-Gaussianity therefore represent one of the 
key tests to detect deviations from the minimal scenario. To 
detect them it is essential to have high sensitivity and large- 
coverage CMB data. 

The Wilkinson Microwave Amsotropy Probe (WMAP) satel- 
lite has produced high resolution all-sky observations of 
the CMB with unprecendented accuracy, which have con- 
firmed many pre dictions of the ACDM concordance cos- 
mological model llSpergel et al. I l2003l : iBennett et al] l2003l : 
ISpergel et aLlbOOel v 

First analyses based on the WMAP first year data and global 
statistical measures did not show significant deviat ions from 
the Gaussianity hypothesis (|Komatsu et al. Il2003h . A num- 
ber of subsequent studies of the first-year WMAP-data 
yielded, however, evidence for non-Gaussian features as well 
as asymmetries in the CMB as measured with the WMAP 
data. 

These results were obtained using a number of qual- 
itatively differ ent statistical tests ranging from (lo- 
cal) spectral (de Oliveira-Costa et al.' 2004; Copi ct al.l 
2004; Erikse n_ct_aL 2004b; Hansen, Banday, & Gorski 2004; 
Landl l2005f) and A^-poin t correl ation function anal- 
ysis llErikseii et all l2004bl. bOOSal ). to Fourier phases 



JChiang et al.l |2003| : 



Coles et al 
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20051). Minkowski - functionals (Park 2004; Erikson ct aLI 
2004al) and wavelet s (tVielva et a l. 2004; Cruz ot al. 2005i: 
McEwen et"al]|2005l ). 

In March 2006 the three-year WMAP results were released 
and several groups revisited the anomalies found in the first 
data release. Considering the fact that the first year WMAP 
data were already of very high quality it was no surprise that 
many of t he anomalies found i n the first year data were re- 
detec t ed (iBrid ges et a l]|2006l: ICopi et al.ll2006l: IJaffe et al.1 



I2OO6I : iLandl I2OO6: Ma rtmez-Gonzalez et al. I bood ) in the 



three year data release albeit none of these tests alone 
showed outstanding significance. Hence there is an ongo- 
ing discussion, whether the identified features are signifi- 
cant enough to prove asymmetry or non-Gaussianity of the 
CMB. Given the tremendous consequences a claim of non- 
Gaussianity of the CMB on cosmology would have, some 
caution in making such statements is warranted. 
Non-Gaussianity cannot be defined in a unique manner, it is 
simply the presence of any higher order spatial correlations. 
Therefore, there cannot be one unique test for all possible 
ways a random field can be non-Gaussian that is prefer- 
able to other tests, or that has superior sensitivity. In or- 
der to perform a thorough analysis, one rather has to apply 
several different, preferably complementary tests. Following 
this reasoning, we present in this paper a novel approach to 
analyse observed CMB-data on a local scale and to test for 
Gaussianity and (statistical) isotropy. Namely we estimate 
the local scaling properties of the maps of tempe rature fluc- 
tuatio n s using weighted scaling indices (see e.g. iRath et al.l 
l|2002l '): lRath fc Schuecker (2003, )). These scaling indices are 
sensitive to the local morphological properties of the field 
of the temperature fluctuations AT{0, (j>) at a given scale r. 
It has been demonstrated that this statistic is highly sen- 
sitive to non-Gaussian signatures in simulated CMB-maps. 
Here, we propose a formalism that can be applied to ob- 
served WMAP-data with its spherical symmetry and test 
for both non-Gaussianties and statistical asymmetries in 



the WMAP three-year data using Monte Carlo simulations 
based on Gaussian fluctuations with a best fit ACDM power 
spectrum and WMAP-like beam and noise properties. 
The paper is organised as follows. In Section 2 we describe 
the WMAP-data as well as the simulations and the prepro- 
cessing steps, which were applied to both the simulated and 
observed data. In Section 3 we review the formalism of the 
weighted scaling indices and explain it with a simple exam- 
ple. As the weighted scaling indices and the wavelets share 
some properties (e.g. locality, scale dependence) we compare 
the two test statistics, show similarities and outline differ- 
ences. We further extend the scaling index formalism to the 
application to spherical WMAP-data. In Section 4 the re- 
sults of the WMAP three-year data analysis are presented. 
We summarize our main results in Section 5 and present our 
conclusions in Section 6. 



2 WMAP DATA, SIMULATIONS AND 
PREPROCESSING 

We use the noise-weighted sum T of the VI, V2, Wl, W2, 
W3 and W4 foreground-cleaned maps. 



T(e,0) = 



(1) 



where j — 5, 6 refer to the two V-band receivers and 
J = 7, . . . , 10 to the four W-bands. a^ j denotes the three- 
year noi se per observat i on fo r the six frequency bands 
given bv lHinshaw et al. I (|2006l ). In contrast to man y of th e 
one-year Gaussian ity studies (e.g. iKomatsu et al. I (|2003l ): 
IVielva et al.1 (|2004h ) the Q bands were no longer taken into 
account, because they exhibit foreground contaminations. 
All these foreground cleaned maps were taken from the pub- 
licly accessible LAMBDA websiteQ- 

In order to study the Gaussianity of the coadded VW- 
WMAP maps, we generated and analysed a set of 1000 
Gaussian Monte Carlo simulations. For this we took the 
set of the best fit Ci's (WMAP only) as published on the 
LAMBDA-site some weeks after the three year data release 
0. For the set of Cj's, random aim's for the Ci's of CMB 
realizations were generated and convolved at each one of 
the WMAP receivers with the appropriate beam transfer 
functions 0. After the transformation from harmonic to real 
space, uncorrelated Gaussian noise realizations were added 
following the number of observations per pixel Nj{9, (j)) and 
the noise dispersion per observation (coj)- We combined all 
the maps following equation (1). Finally, both the coad- 
ded WMAP-data and the simulations were degraded to 
Nside ~ 256, the KpO mask was applied, and the residual 
monopole and dipoles were fitted and subtracted. All these 
simulation and prepr ocessing steps were performed with the 
HEALPIX-software iGdrekTeiaDllooi). 



^ http://lambda.gsfc.nasa.gov 

^ http: / /lambda. gsfc.nasa.gov/product / map / current / params/lcdm_wmap.cfm 
http: / /lambda. gsfc . nasa.gov / product /map/dr2 /xfer_funcs_got . cfm 
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3 WEIGHTED SCALING INDICES 



The weighted scaling ind ex method (SIM) (|Rath et al.ll2002l : 
iRath fc Schuecke j l2003h offers one possibihty to estimate 
the local scaling properties of a set of points, which is gen- 
erally repesented in a d— dimensional space. The SIM has 
found applications in (astrophysical) time series analysis of 
AGNs, where the observed light curves were represented as 
point distributions by embedding the time series in a higher 
dimensional a rtificial phase space using delay coordinates 
jGliozzi. et al.,,2002. : .Gliozzi, Papadakis, fc Rath.200^ ') . 
Scaling indices have also often successfull y been used for 
structure analysis in 2D a.nd 3D image data |jamitzkv et alJ 
I2OOII : iMonetti et "alll2003l : iMueller et al.ll2004l '). For this the 
image data are represented as point distributions by com- 
prising the spatial and intensity information of each pixel. 
For two-dimensional images one thus obtains a set of three- 
dimensional vectors pi = {xi, yi, I{xi, yi)),i = 1, . . . , Npi^eis- 
On the basis of the representation of images as point dis- 
tributions the weighted scaling indices are calculated as fol- 
lows: 

3.1 General Formalism 

Consider a set of N points P = {j>i\, i — 1, . . . , N . For each 
point the local weighted cumulative point distribution p is 
calculated. In general form this can be written as 



Pifi 



(2) 



where Sr{») denotes a kernel function depending on the scale 
parameter r and d{») a distance measure. 
The weighted scaling indices a{pi,r) are obtained by calcu- 
lating the logarithmic derivative of p(pi,r) with respect to 
r, 

a(pl..)- "y-^^ --^p(pl,r). (3) 
O log r p or 

In principle, any differentiable kernel function and any dis- 
tance measure can be used for calculating a. In the following 
we use the euclidean norm as distance measure and a set of 
gaussian shaping functions. So the expression for p simplifies 
to 



pifi 



N 



di 



(4) 



The exponent q controls the weighting of the points accord- 
ing to their distance to the point for which a is calculated. 
The higher q is the more steplike becomes the weighting 
function resembling more and more the Heaviside-function, 
which is used for the calculation of the unweighted scaling 
indices. Another interesting choice of q is given by g = 2. 
In this case the kernel function is the well-known Gaussian 
exponential function. 

Throughout this study we calculate a for the case q ~ 2. In- 
serting the expression (4) in the definition for the weighted 
scaling indices in (3) yields after some algebra the following 
analytical expression for a: 



a{Pi,r) = 



.7 = 1 



(5) 



Structural components of a point set are characterized by 
the calculated value of a. For example, points in a point-like 
structure have a ~ and pixels forming line-like structures 
have a ~ 1. Area- like structures are characterized by a « 2 
of the pixels belonging to them. A uniform distribution of 
points yields a « d which is equal to the dimension of the 
configuration space. The scaling indices for the point set 
form the frequency distribution N{a) 



N{a)da = #(q £ [q, a + da[) 

or equivalently the probability distribution 

P{a)da = Prob{a £ [q, a -I- da[) 



(6) 



(7) 



The P(a)-representation of a point set can be regarded as 
a structural decomposition of the data where the points are 
differentiated according to the local morphological features 
of the structure elements to which they belong to. 



3.2 Scaling Indices and Wavelets 

Some test statistics, with which non-Gaussian signatures 
have been detected in the first and three-year WMAP- 
data . are based on a wavelet analysis of the CMB data, 
e.g. (IVielva et al ]|2004l : lMcEwen et al.ll200 5l). Beside the fact 
that wavelet-based tests yield one of the highest significances 
for non-Gaussianities in the WMAP-data, wavelets and scal- 
ing indices can also be considered to be similar to each other. 
Both measures perform a local analysis of e.g. image data 
and they are calculated for different scales, which yields 
information about characteristic sizes of detected features. 
On the other hand, there are obvious differences. While the 
scaling indices are a nonlinear and non-bijective local filter, 
the wavelet transformation is a linear operation, namely a 
scale- and space-dependent filtering of the image data with 
a wavelet function 



dxdyl[x,y)^{x,y-b^,hy,i 



(8) 



For the special case of spherical mexican hat wavelets 
(SMHW), the 'mother wavelet' ilj{d,r) is defined by: 



vE'(d,r) 



with d 



*(d,6, ^0,by =0,r) 



(27r)i/2 



(9) 



^. It has been shown 

(jMartmez-Gonzalez et al. II2002I ') that SMHW yield superior 
results in detecting non-Gaussian signatures compared to 
other filter functions, e.g. Haar wavelets. 
To consider the performance of the scaling indices in the 
context of known test statistics for non-Gaussianities we 
calculate and compare scaling indices and SMHW for two 
examples. 

In the first simple example we generated a synthetic image 
which consists of a (white Gaussian) noisy background, 
where the noise level was chosen to be twice as high in 
the lower third of the image than in the rest of the image. 
We interspersed three single- coloured lines and six disc-like 
elements in the noisy background. The intensities of the 
structural elements ranged from Oct to ±la with respect to 
the standard deviations of the different noisy backgrounds 
(see Figure [T]). 
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Figure 1. Left to right: Syntliotic test image, response image of 
the scaUng indices and response image of the wavelet coefficients. 
The response images were normalised to the range of values of 
the scaling indices and wavelets respectively. 
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Figure 2. Probability distribution of the the scaling indices 
(above) and wavelet coefficients (below). The gray lines indicate 
the distributions for the pure noise image without the line- and 
disc-like structural elements. 



We calculated the wavelets and scaling indices for a cor- 
responding and well-suited scale. Figure [T] shows the grey- 
value coded response of the wavelet and scaling index fil- 
ter. One can see that the wavelets can only detect those 
discs, which deviate most from the mean value of the noise. 
The lines are only hardly detected in the less noisy back- 
ground. The scaling indices, however, can detect all inter- 
spresed structural elements in the image, where the (inner 
part of the) disc-like structures have lower a- values than the 
three lines. The observed differences in the images of the 
filter response find their refiections in the probability distri- 
butions P{ot) and P(w) for the scaling indices and wavelet 
coefficients (Figure |2]). 




Figure 3. Simulated non-Gaussian CMB map with superimposed 
white additive noise (lower left) and three surrogate realisations 
with the same power spectrum and amplitude distribution. 
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Figure 4. Significances of the moments of the distribution of scal- 
ing indices (black) and wavelet coefficients (gray) for the image 
in fig. [3] The lines and + denote the mean, the dotteds line and * 
the standard deviation, dashed lines the skewness and dashdotted 
lines and o the kurtosis. 



While the spectrum of scaling indices shows a clear sec- 
ond peak at a « 2.1 and a larger tail towards lower a- values, 
no visible deviations can be seen in the in the distribution 
of the wavelet coefficients when compared with those for a 
pure noisy image. With this example it becomes obvious 
that the wavelets are more sensitive to structures, which are 
associated with intensity variations of significant magnitude 
with respect to image noise. The scaling indices can, how- 
ever, detect also structural features, which do not manifest 
themselves with significantly higher (lower) intensity values 
but as intrinsic structural variations within the noise level. 
In the second more realistic example we consider a realisa- 
tion of CMB anisotropies due to the Kayser-Stebbin effect 
from cosmic strings on which white additive Gaussian noise 
was superimposed. The noise level is chosen with rms signal- 
to- noise-ratio (S/N) of 0.5. For this image 20 surrogate maps 
were generated, which have the same power spectrum and 
amplitude distribution as the original simulation (Figure IS)) • 

For more details about th e simulation and its surrogates 
see iRath fc Schueckerl l|2003l ). We compare tests for non- 
Gaussianities based on wavelets and scaling indices. In both 
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cases we use the moments of the distributions of wavelet 
coefficient and scahng indices as test statistics, namely the 
mean (< a >), standard deviation (tJa), skewness (7a) and 
kurtosis (k^), which are defined by 



f 

1=1 



N 



1\ 



1/2 



< a > 



i = l 



7c« 



AT , 

_L \ " ( 

N 2^ 

i = l 
N , 

N ^ 

i = l 



< a > 



< Q > 



(10) 



(11) 



(12) 



(13) 



as test statistics. For each statistic the degree of de- 
tected non-Gaussianity is measures in terms of significance 
5, 



M- <M> 



Cm 



(14) 



where M is < a >,aa,"/ci and Ka respectively. The mean 
and standard deviation of the measure M are derived from 
the 20 surrogate reliasations. Figure 2] summarizes the re- 
sults for the scaling indices and SMHW. Both classes of fil- 
ters can well detect the non-Gaussianities even at this high 
noise level. The highest significances occur at similar scales. 
However, the best results were obtained for different mo- 
ments. While for the scaling indices the standard devaition 
performed best, we found for the wavelets the best discrim- 
ination with the kurtosis. The second best result for the 
scaling indices was obtained with the mean. Based on these 
findings we restrict ourselves to the analysis of the mean and 
standard deviation of the distribution of scaling indices for 
the WMAP data. In summary we found in this example sim- 
ilar discrimination results for the two classes of filters, but 
for different moments. The scaling indices gave - in general 
- better results for the lower moments of the distribution 
than the wavelets. Since the lower moments are less sen- 
sitive to outliers of the distribution, one can arguably say 
that an analysis of the first moments of the scaling index 
distribution is statistically more stable and thus preferable. 



3.3 Weighted Scaling Indices for WMAP-Data 

To apply the SIM to the WMAP-data one has to find a 
proper representation of the spherical data in a suitable em- 
bdding space. Our aim in this work was to stay as close as 
possible to the original data, i.e. to maintain the spherical 
symmetric character of the data set and thus to omit any 
projection onto a flat space. On the other hand, we strove 
to find a proper three-dimensional embedding of the data 
according to the three free parameters of each pixel on the 
unit sphere, namely their angles and temperature. For this, 
we chose the following approach: 

Consider a temperature map I (6, (j)), where each pixel is as- 
signed with a temperature value I{6,(f)). The position of 
the pixel on the (unit) sphere is given by the two angles 







7^1 
































Figure 5. Left: WMAP-data represented as a three-dimensional 
point distribution. Right: x, z-projection of all points with \y\ < 
0.05. Above: a = 0.075, below: a = 0.225. The black circles indi- 
cate the scaling ranges r = 0.075 and r = 0.225 respectively. A 
good sensitivity to the temperature fluctuations at a given scale 
r is obtained for r = a. 



9 and <j). In addition, it contains the temperature informa- 
tion I(9,(f>). One possible three-dimensional representation 
of the WMAP-data, in which both the spatial and temper- 
ature information of each pixel is simultaneously taken into 
account, is given by 



X = {R + dR) cos{(j>) sm{e) 
y = {R + dR) sm{(t>) sm{e) 
z = {R + dR)cos{e) 

with 

dR = a{T{e,(t))- < T >)/(JT . 



(15) 
(16) 
(17) 

(18) 



By introducing the term dR the temperature 
anisotropies are transformed to variations in the radial di- 
rection around the sphere. The normalisation ensures that 
dR has zero mean and a standard deviation of a. The pa- 
rameters R and a are free (scale) parameters, which control 
the size of the temperature-induced radial jitter relative to 
its spatial extent. These two parameters R and a, together 
with the scaling range parameter r for the calculation of the 
scaling indices have to be properly set so that the scaling 
indices are sensitive to the temperature fluctuations in the 
chosen (spatial) scaling range (see Figure [5]). 

To achieve this, we coupled the parameter a, which con- 
trols the width of the radial jitter with the scaling range r 
and set for each scaling range a = r. The size of the the 
sphere R has to be chosen large enough with respect to the 
radial jitter dR, so that the scaling indices remain a local 
measure and their calculation are not affected by pixels on 
the opposite side of the sphere with a large dR. Hence, we 
deliberately set R — 2 for aU calculations presented below. 
We performed some tests for the WMAP data and a subset 
of simulations using different values of R, e.g. R = 1.5, and 
found only marginal differences in the results. 
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Figure [5] shows the representation of the WMAP data as 
a point set P = {pi = {xi,yi, Zi)},i = l,...,NpM and 
two-dimensional projections of all data with \y\ < 0.05. For 
each point pi a set of weighted scahng indices a for ten 
different values of r, r = 0.025, 0.05, 0.75, . . . , 0.25 is calcu- 
lated, which (roughly) corresponds an angular resolution 
of ^ = 1.4°, 2. 8°,..., 14.3°. 

Thus, the weighted scaling indices can also be interpreted 
as a filter response of a local nonlinear filter acting on the 
spherical CMB data. Global quantities of the Q-distribution 
as well as the probability density P{a) are used to derive 
statistics quantifying the structures in the WMAP-data and 
respective simulations. 



4 RESULTS 

In Figure IH] the preprocessed foreground-cleaned coadded 
map of the WMAP data and the a- response for two differ- 
ent radii r is shown. It is obvious that for these larger scales 
the pixels in the vicinity of the KpO-mask are affected by 
boundaries effects, which lead to systematically lower values 
for a. This effect, which is only due to the lack of data points 
within the KpO mask, is the same for both the WMAP-data 
and the simulations, so that global measures based on the 
scaling indices, e.g. the moments of the probability distribu- 
tion, are systematically affected in the same way. However, 
some smaller local effects may be diluted by the boundary 
effects. In order to have a cleaner map we also considered 
only scaling indices in a more conservative mask, where only 
pixels with |6| > 30 (6: Galactic latitude) were selected (Fig- 
ure[5]). Note that for the calculation of the scaling indices all 
pixels outside the KpO-mask were taken into account. In the 
following we always explicitly specify to which chosen mask 
the presented results refer. 

The probability densities P{a) of the scaling indices 
for one selected scale {r — 0.175) and both masks are dis- 
played in the Figures [7] and [8] for the WMAP-data and a 
subset of 20 simulations. One immediately realizes that for 
both masks the probability density for the WMAP data is 
shifted towards higher values, which indicates that the un- 
derlying temperature fluctuations for the observed data re- 
semble more 'unstructured', i.e. purely random and uniform 
fluctuations than the simulations. This effect is more pro- 
nounced in the northern hemisphere than in the southern. 
Furthermore, the distributions are broader for the simula- 
tions than for the WMAP-data, indicating that the simula- 
tions exhibit a larger structural variability than the observed 
data. 

These effects can more rigorously be quantified by calcu- 
lating the mean and standard deviation for the distribution 
of scaling indices as calculated for different scaling ranges. 
Figures l9l and [TOl show these results. For scales larger than 
r = 0.1 the mean of the scaling indices is always systemat- 
ically higher for the WMAP than for the simulations. The 
effect is much more pronounced in the northern hemisphere. 
For the standard deviation we observe for the same scales 
significantly lower values for WMAP in the northern hemi- 
sphere and slightly higher ones for the southern sky. For the 
full sky these two effects cancel each other so that the ob- 
served deviations to lower values are no longer so significant. 
Beside the mean and standard deviation we additionally 






Figure 6. From top to bottom: Co-added WMAP-map with KpO- 
mask, color-coded a-response for r = 0.175 and r = 0.225 and 
KpO- and extended mask. 



considered a combination of these two test statistics, namely 
a diagonal x^-statistic 

(TMi 

where M\ =< a > and M2 = (J a- This statistic is com- 
puted for both the simulations and t he observed moments . 
Note that we follow the reasoning of lEriksen et all (|2004ah 
and choose a diagonal x^-statistics, because also in our case 
the moments are highly correlated leading to high values 
in the off-diagonal elements of the cross-correlation matrix. 
But if the chosen model is a proper description of the data, 
any combination of measures should yield statistically the 
same values for the observations and the simulations. The 
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Figure 7. Probability density P{a) for the co-added WMAP 
data and 20 simulations for r = 0.175 and the KpO-mask. 
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Figure 9. Scaling index mean statistics as a function of the scal- 
ing range r for the full sky (above), northern (middle) and south- 
ern hemisphere (below) using the KpO-mask. -|- denote the values 
for WMAP, the dark gray and light gray areas indicate the Icr 
and 3cr regions around the mean value (black line) as derived from 
the simulations. 



Figure 8. Same as Figure [7] but for the extended mask. 
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Figure 10. Same as fig. |9] but for the standard deviation a^- 
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Figure 11. Significances of tlie (combined) moments of the a- 
distribution of the WMAP-data with the KpO-mask as a function 
of the scaling range r. + denotes the mean, * the standard de- 
viation and the boxes x^-combination of the mean and standard 
deviation. 



significances of the deviations of the WMAP-data from the 
simulations 

M- <M> 



S = 



0"M 



(20) 



(M =< a >,a„ and y^) are shown in Figures [TT] and [T^ for 
the KpO-mask and the extended mask. 

The mean < M > and standard deviation cjm of the 
three measures M were derived from the set of simulations. 
As a pure frequentist approach, we also considered the con- 
fidence or significance level of the null hypothesis that the 
observation belongs to the Gaussian Monte Carlo ensemble 
and consider the fraction of simulations p, which have higher 
(lower) values of M than the moments for the observation. 
In the Tables [TJ [5] and [3] the significances S as well as the 
probabilities p are listed. 

Analysing the significances and the confidence levels of 



Figure 12. Same as Fig. Illl but for the extended mask. 



the (combined) moments of the distribution of scaling in- 
dices as a function of the scaling range the results read as 
follows: For the full sky we obtain for the mean value < a > 
significances ranging from 2.0 to 2.2 for the different scales. 
The confidence levels for the detection of non-Gaussianties 
are, however, very high and do not fall below 99% for any 
scale. A closer look at the distribution of < a > for the sim- 
ulations reveals that this distribution is not Gaussian but we 
find a number of outliers with very high values for < a >, 
which leads to high (T<a> thus to low significances. So al- 
though the significances for detecting non-Gaussianities are 
not very high, the confidence levels, for which no implicit as- 
sumptions about the distribution function of the test statis- 
tics is made, strongly indicate the presence of non-Gaussian 
features in the observed WMAP-data. Even higher values for 
both the significances and the confidence levels are found, if 
one only considers the northern hemisphere. In this case the 
sigificances range from 1.9 for smaller scales to 2.8 for the 
largest scale. For scales larger than r = 0.15 none of the sim- 
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Table 1. < a > with KpO mask 



Table 4. , all scales 



Scaling Range 


Full Sky 


Northern Sky 


Southern Sky 


) 




(a/ /o) 


(a/ /o) 


0.025 


2.2/99.5 


2.1/99.4 


2.3/99.8 


0.050 


2.1/99.8 


2.0/99.3 


2.2/99.7 


0.075 


2.0/99.5 


1.9/99.3 


2.0/99.4 


0.100 


2.0/99.4 


2.0/99.7 


1.9/98.8 


0.125 


2.0/99.5 


2.1/99.8 


1.9/98.6 


0.150 


2.1/99.6 


2.4/99.9 


1.8/98.4 


0.175 


2.2/99.7 


2.5/>99.9 


1.7/97.4 


0.200 


2.2/99.7 


2.7/>99.9 


1.6/96.0 


0.225 


2.2/99.7 


2.7/>99.9 


1.5/94.0 


0.250 


2.1/99.4 


2.8/>99.9 


1.3/91.6 


Table 2. a a with KpO mask 






Scaling Range 


Full Sky 


NortrliGrii Sky 


SoutliGrn Sky 


(r) 


iS/%) 


(S/%) 


(S/%) 


0.025 


2.0/98.9 


1.9/98.7 


1.9/98.9 


0.050 


2.2/99.7 


2.0/98.9 


2.1/99.1 


0.075 


0.4/64.8 


0.2/60.3 


0.4/63.8 


0.100 


1.6/94.9 


1.5/94.2 


0.9/81.6 


0.125 


2.2/99.0 


2.6/99.8 


1.0/84.6 


0.150 


2.3/99.0 


2.9/99.9 


0.8/79.5 


0.175 


2.1/98.4 


2.9/99.9 


0.6/71.6 


0.200 


1.9/97.5 


2.8/99.9 


0.5/66.8 


0.225 


1.7/95.8 


2.6/99.7 


0.4/64.3 


0.250 


1.5/93.1 


2.4/99.7 


0.2/58.3 



ulations was found to have a higher values for < a > than 
the observation, which represents a quite unambigous detec- 
tion of non-Gaussianties in the northern hemisphere. For the 
southern hemisphere, however, both the significances and 
the confidence levels for the smaller radii are slighly higher 
than for the northern sky but continously decrease for higher 
radii r. For the standard deviation we find slighly different 
results. For the smallest scales (r — 0.025) is signifi- 
cantly larger for WMAP than for the simulations. In a tran- 
sition regime r « 0.075 the standard deviation is practically 
the same for the observation and the Monte Carlo sample. 
For larger scales we observe higher standard deviations for 
the simulations. This effect is much more pronounced in the 
northern hemisphere giving rise to significances up to 2.9 
and very high confidence levels for intermediate scales. For 
the largest scales the differences for a between simulations 
and oberservation diminishes. For the southern hemisphere 
the width of the distributions become more and more simi- 



Table 3. with KpO mask 
Scaling Range Full Sky 



Northern Sky Southern Sky 



(r) 


is/%) 


(■S/%) 


is/%) 


0.025 


2.0/96.7 


1.8/95.9 


2.1/96.9 


0.050 


2.3/97.5 


2.0/96.1 


2.4/97.7 


0.075 


0.9/88.4 


0.7/85.9 


1.0/88.4 


0.100 


1.8/95.4 


1.7/94.9 


1.1/89.3 


0.125 


2.7/97.8 


3.6/98.9 


1.0/89.2 


0.150 


2.9/98.2 


4.8/99.3 


0.8/85.7 


0.175 


2.7/97.9 


5.2/99.4 


0.5/80.3 


0.200 


2.5/97.1 


5.1/99.3 


0.3/76.4 


0.225 


2.1/96.6 


4.8/99.3 


0.1/70.1 


0.250 


1.8/94.7 


4.5/99.0 


0.1/63.0 



Full Sky 

iS/%) 



Northern Sky 

(S/%) 



Southern Sky 

(S/%) 



2.1/96.9 
2.4/96.5 
2.4/97.3 



2.7/97.7 
4.4/99.5 
3.7/98.9 



1.5/94.2 
0.2/70.0 
1.1/91.6 



lar so that no signatures for deviations from Gaussianity are 
identified using aa- The behaviour of the x^- statistics as a 
function of the scale parameters r can - as expected - be 
regarded as a superposition of the two underlying statistics 
< Of > and Ga- Only the significances are highly increased 
leading to a 5a detection of non- Gaussianity at scales of 
r = 0.175 in the northern sky. Performing the same analy- 
ses for the extended mask (Figure [T^J yields essentially the 
same results. Only marginal variations for the significances 
and confidence levels (which are not explicitly shown in this 
paper) are found in this case. Thus, these results are quite 
stable with respect to differences for the chosen mask. 
Some readers might argue that the selection of the mo- 
ments and especially of the scales, where the significances are 
largest, represents an a posteriori choice taken after looking 
at the data. However, the choice of the moments was moti- 
vated by the results obtained with simulations (see Section 
3) and the investigation of the data on different length scales 
can be regarded as a fairly conventional and unbiased ap- 
proach. If a random field is Gaussian, it must be Gaussian 
on all scales. If deviations are detected at some scales, one 
can already infer non-Gaussianity. Nevertheless, we also cal- 
culated diagonal x'^-statistics, where we considered only one 
(mean or standard deviation) or both measures, and summed 
over all length scales , 



JV,, 

i = l 
2 \ ^ 



A/i(rO- < Afi(rO > 



=1 



Mo 



0"Mi(ri) 



< Ma 



0"Mi(ri) 



Nr 2 

EE 



(21) 



(22) 



(23) 



with Afi(ri) =< aifi) >, M2(ri) = (Ja{ri) and A^,- being 
the number of considered length scales (here: Nr = 10). The 
results are summarised in Tabled] 

Also for this obviously a priori test statistics, where 
some unimportant scales contribute to the final value of x^ 
and may dilute the result, we find significant signatures for 
non-Gaussianities in the northern sky. 

Beside the previous analyses based on global measures de- 
rived from the P(ot)- distribution, we also investigate the 
spectrum of scaling indices on a differential level. Therefore 
we consider the bin- wise significances S{oii), 



P(qO- < P(qO > 



(24) 



for the probability densities P{a) of selected scales (see Fig- 
ures [T3] and [TH) . 

Note that we omitted in this case the absolute value in 
the definition of the significance in order to show in which 
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Figure 13. Bin-wise significances (with sign) as derived from tlie 
pribability distribution P{a) for the WMAP data and KpO-mask 
for r = 0.175 (above) and r = 0.225 (below). 
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Figure 15. Marked pixels with a g [2.0,2.3] (yellow), a 6 
[2.425,2.475] (blue) and a G [2.60,2.65] (red) for r = 0.175 
(above) and r = 0.225 (below). 

direction with respect to the simulations the P{a) of the 
WMAP data deviates. Although the spectra for the KpO- 
mask (Figure I13|l are diluted by edge effects for smaller a, 
one can clearly detect the sharp transition from systemati- 
cally small to high values in the significances, which is due 
to the shift of the whole spectrum towards higher a-values 
for the WMAP-data. This shift is so pronounced that for 
some bins of the spectrum even the 3a- level for the de- 
viation is nearly reached. If we consider the same spectra 
for the extended mask, where almost all disturbing edge ef- 
fects are removed, we see (Figure [14]) - beside the transition 
from negative to positive values due to the global shift - a 
new, highly significant feature at a « 2.1 in the southern 
sky emerging, when the scaling range is increased. There is 
an excess of pixels with these low a— values in the WMAP- 
data as compared to the simulated maps. To elucidate the 
origin of these pixels we color-coded the pixels according 
their value of scaling indices (Figure [15)). 
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Figure 14. Same as fig. I13l but for the extended mask. 
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Figure 16. Significancs of the (rotated) upper hemisphere for 
the SIM based statistics for r = 0.175. From top to bottom: 
mean < a >, standard deviation cfol and x^-combination of mean 
and standard deviation. The highest significances were obtained 
for {e,ip) = (27°, 35°) for < « >, {e,<f>) = (39°, 45°) for (Tq and 
{0,<f>) = (30°, 41°) for ■ The maximal asymmetry between the 
rotated northern and southern hemisphere were found for the 
rotation angles {9,(f)) = (27°, 35°) for < a >, (6»,0) = (30°, 41°) 
for o-Q and (9,0) = (30°, 41°) for x^. 





Figure 17. Same as Fig. \W\ but for r = 0.225. For this scal- 
ing range, the highest significances were obtained for {9, 0) = 
(21°, 148°) for < a >), (e,4>) = (27°, 65°) for CTq and (9,4>) = 
(30°, 41°) for • The maximal asymmetries between the rotated 
northern and southern hemisphere were found for the rotation 
angles {e,(f>) = (21°, 148°) for < a >, {e,<j}) = (27°, 55°) for era 
and (e,0) = (30°, 41°) for x^- 
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It becomes immediately obvious that the pixels with 
the small a- values identified in the southern sky form two to 
three clusters (depending on the scale r), where the largest 
one corresponds to the cold spot - the well-known anomaly 
first detected bv lVielva et al.l l|2004 ) by means of a wavelets 
analysis. The pixels belonging to the other two a-intervals, 
for which the largest deviations in the bulge of the distribu- 
tion are found, cannot be associated with special localised 
features. These pixels are rather distributed all over the sky, 
which indicates that the shift of the spectrum of scaling in- 
dices represents a global intrinsic effect. 
The distinction between northern and southern hemisphere, 
which we made during the previous analysis, is natural be- 
ing triggered by excluding the foreground-contamined area 
of the galactic plane. On the other hand, this choice is arbi- 
trary, because no symmetry axis is preferable for a presum- 
ably isotropic CMB. To test for asymmetries in WMAP- 
data we consider rotated hemispheres, calculate the global 
measures < a >,aa and for 3072 rotation angles and 
compare the results for the (rotated) northern and southern 
hemispheres. In Figures [16] and [IT] we show the significances 
S{< a >), S{aa) and 5(x^) of the northern hemisphere as 
determined in a reference frame where the north pole pierces 
the center of the color-coded pixel. All three measures yield 
systematically higher significances for rotations pointing to 
northern directions relative to the galactic coordinate sys- 
tem. The rotation angles, for which the highest significances 
of the northern hemisphere are obtained, are listed in the 
captions of the respective figures. They are very similar but 
do not always coincide with the rotations, for which the max- 
imal asymmetries, as measured by the difference of the sig- 
nificances of the rotated northern and southern hemisphere 
are found. It is worth noticing that the direction of the most 
pronounced non-Gaussianities and asymmetries differs from 
the dipole direction {9,4i) = (42°, 264°) and the so-called 
axis of evil {6,<j)) ~ (30°, 260°), which is very close to the 
dipole direction. 



5 SUMMARY 

We performed a scaling index analysis of the WMAP three- 
year data. Specifically, we analysed the foreground-cleaned 
co-added maps of the V-and W-band. We found highly sig- 
nificant signatures of both non-Gaussianities and asymme- 
tries in the WMAP three-year data. Our main findings can 
be summarized and interpreted as follows: 
In the northern hemisphere the spectrum of scaling indices 
is systematically broader and shifted towards higher values 
yielding highly significant deviations of the mean and stan- 
dard deviation of the distribution. This effect can naturally 
be interpreted as too few structure and structural variations 
in the CMB-fiuctuations as measured by WMAP compared 
to the predicted ones within the concordance model. The 
highest global signatures for non-Gaussianities and asym- 
metries between the northern and southern hemisphere were 
found for rotated coordinate systems, where the significances 
for the detection of signatures for non-Gaussianities range 
from 2.6a up to 7.4a. These findings are q uite consistent 
with p re vious result s for the first year data (|Eriksen et al.l 
l2004al lbl: iPark II2004I ) , where very similar features were iden- 
tified - even though with a smaller significance level - using 



the Minkowski functionals, the power spectrum and A'^-point 
correlation functions. 

In the southern hemisphere the global properties of the P{a) 
distribution for the WMAP-data are more consistent with 
the simulations than in the northern sky. We find, however, 
highly significant localisable features of non-Gaussianity, for 
which the largest one can be as sociated with a col d spot. 
An anomalous area detected by IVielva et al] (|2004!) in the 
WMAP fir st year data and confirmed in the WMAP three 
year data (|Cruz et al.ll2006l '). 



6 CONCLUSIONS 

In conclusion, we demonstrated for the first time the feasi- 
bility to adapt and apply the scaling index method as an 
estimator of local scaling properties of a point set to spher- 
ical symmetric data. 

The results obtained with the scaling indices give further 
strong evidence that also the coadded WMAP three year 
data do indeed contain unusual features, which are not in 
agreement with the hypotheses of Gaussianity and isotropy 
predicted in the standard inflationary scenario. 
A quite remarkable result of our study is the fact that two 
previously known anomalies, namely the lack of structure in 
the northern sky and the cold spot in the southern hemi- 
sphere could be reidentified using a completely different test 
statistic. This increases the evidence that these anomalies 
are of true physical origin and not related to one single test 
statistic. 

It is also worth noticing that the scaling indices could de- 
tect both of these anomalies simultaneously, whereas the 
wavelets were only sensitive to localised spot-like struc- 
tures and the A''-point correlation functions and Minkowski- 
functionals only detected the non-Gaussianities and asym- 
metries on large scales. 

The main task for future studies is to elucidate the possible 
sources of these anomalies, whether they are due to sys- 
tematics or foreground effects or truely represent intrinsic 
CMB-fluctuations due to some exotic physics. Future more 
detailed studies will investigate the possible origins of the 
non- Gaussian signatures by separately analysing the V and 
W bands and by comparing the three year WMAP-data with 
the first, second and third year data. All these tests may 
help to find hints about the origin of the non-Gaussian sig- 
natures. 

In either case we could demonstrate that the scaling index 
method provides very flexible and highly sensitive statistics 
for e.g. the identification of asymmetries and non-Gaussian 
signatures in the WMAP-data thus representing an impor- 
tant novel statistical tool for CMB analyses. 
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